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The Causes of Halo Shape Changes Induced by CooUng Baryons: 

Disks Versus Substructures 

Victor P. Debattista^'^, Ben Moore^, Thomas Quinn^, Stelios Kazantzidis^, Ryan Maas^, 
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ABSTRACT 

Cold dark matter cosmogony predicts triaxial dark matter halos, whereas ob- 
servations find quite round halos. This is most likely due to the condensation 
of baryons leading to rounder halos. We examine the halo phase space distribu- 
tion basis for such shape changes. Triaxial halos are supported by box orbits, 
which pass arbitrarily close to the density center. The decrease in triaxiality 
caused by baryons is thought to be due to the scattering of these orbits. We 
test this hypothesis with simulations of disks grown inside triaxial halos. After 
the disks are grown we check whether the phase space structure has changed 
by evaporating the disks and comparing the initial and final states. While the 
halos are substantially rounder when the disk is at full mass, their final shape 
after the disk is evaporated is not much different from the initial. Likewise, the 
halo becomes (more) radially anisotropic when the disk is grown, but the final 
anisotropy is consistent with the initial. Only if the baryons are unreasonably 
compact or massive does the halo change irreversibly. We show that the char- 
acter of individual orbits is not generally changed by the growing mass. Thus 
the central condensation of baryons does not destroy enough box orbits to cause 
the shape change. Rather, box orbits merely become rounder along with the 
global potential. However, if angular momentum is transferred to the halo, ei- 
ther via satellites or via bars, a large irreversible change in the halo distribution 
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occurs. The ability of satellites to alter the phase space distribution of the halo 
is of particular concern to galaxy formation simulations since halo triaxiality can 
profoundly influence the evolution of disks. 

Subject headings: galaxies: evolution — galaxies: formation — galaxies: halos 
— dark matter 



1. Introduction 



The dark matter halos that form via hierarchical growth in the cold dark matter (CDM) 
cosmologies are generally triaxial with mean axial ratios 6/a ~ 0.6 a nd c/a ~ 0.4 , where 



c < b < a are the short, intermediate, and long axes, respectivelv (IBardeen et al. 
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sensus from these studies is that dark matter halos are rounder than those predicted by 
collisionless CDM simulations. But this need not be in disagreement with CDM since the 
condensat i on of baryons to the centers of halos has be en shown to lead t o rou nder halos 
(iDubinskil Il994j : iKazantzidis et al.l 120041 ). For example, iKazantzidis et al.l (120041 ) find that 
the principal axis ratios increase by ~ 0.2 — 0.4 in the inner regions (although triaxiality is 
not completely erased) extending to almost the virial radius. 



Slowly rotating triaxial structures can be supported by centrophilic box orbits (jSchwarzschild 



1979 
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). These scattering events lead to a large number of orbits becoming chaotic. Chaos 



by itself, however, need not be a fundamental limit to forming long-lived triaxial struc- 
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tures: using orbit superposition, iPoon fc MerrittI ((2002) were able to construct long-lived 
triaxial models of nuclei even in t he presence of a large frac tion (> 50%) of chaotic orbits. 
If axisymmetrization does occur, I Gerhard fc Binneyl (119851 ) predict that it is largely con- 
fined to the center and o ccurs gradually. The iV-body simulations of a cored system by 
Merritt fc QuinlanI (119981 ) instead found that the axisymmetrization extends to the entire 
system and occurs on a crossing ti me for black holes of mass ^ 2 % of the galaxy's mass. 



When instead the system is cuspy, iHolley-Bockelmann et al.l (120021 ) found that black holes 
do not lead to a global axisymmetrization of the system. Triaxial structures in disks {i.e. 
bars) can also be destroyed by central mass concentrations (CMCs). The main mechanism 
is ag ain scattering by th e CMC. Although the main bar-supporting orbit family, the Xi or- 
bits (jContopoulod Il980l ) . is a centrophobic loop family, stars librating about the closed xi 
orbits can still get close to the center and then be scattered by a CMC. Simulations have 
shown that the required mass for a soft CMC {i.e. one with a scale of a few 100 pc) is an 
unrealistically large ~ 20% of the disk mass, w hile the mass required of a hard CMC (few 
parsecs or less scale) is ~ 5% of the disk mass (IShen fc SellwoodI l2004l : lAthanassoula et al. 



20051 : iDebattista et al.ll2006l ). which is much larger than typical supermassive black holes. 



Likewise, it has often been assumed that the loss of halo triaxiality when baryons cool 
inside halos is partly or mostly due to the destruction of box orbits, which pass arbitrarily 
close to the center after a sufficiently long time. The fate of box orbits in the presence 
of disks is of interest for various reasons beside the shape of the halo. Box orbits play an 
importan t role in speeding up the mergers of supermassive black holes at the centers of 
galaxies (IMerritt fc Poonll2004l ). Moreover, box orbits lead to radial anisotropy, whereas 
the destruction of box orbits results in tangential anisotropy. This in turn affects the event 
rate and energies of dark matter detection experiments involving both direc t scattering and 



indirect annihilation from capture by the Sun or the Earth (see the review by lJungman et al. 



19961). 



In order to help shed light on these issues, we test whether box orbit scattering is 
responsible for triaxial halos becoming rounder when baryons cool inside them. We do 
this via simulations in which we first grow and then evaporate disks inside triaxial halos. 
Such evaporation, while obviously unphysical, allows us to directly assess the impact of 
disks on halos by comparing the initial and final states. After the disks are grown, we find 
that the halos become substantially rounder and their kinematics radially anisotropic. But 
comparing their initial and final shapes when the disk mass is zero in both cases, we find 
that the changes are largely reversible. The destruction of box orbits being irreversible, halos 
should not recover their initial states if this is the main cause of the shape change. We also 
show that if angular momentum is transferred to the halo (via bars or satellites), then the 
irreversible changes are substantially larger. 
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Section [2] of this paper discusses the A^-body methods used in this study. Section [3] 
presents the resuhs of simulations with growing rigid central massive objects. In §1] we 
present simulations in which angular momentum is transferred to the halo either by a live 
bar or by satellites. Section O presents a preliminary analysis of the orbital evolution for a 
subsample of the simulations. Our conclusions are presented in ^ 



2. Numerical Methods 

The basis of this work is that box orbit destruction is an irreversible process. Rather 
than following all orbits, we adiabatically grow and then evaporate a disk to show that 
the distribution function of a halo is not substantially changed despite the fact that the 
halo appears very different when the disk is at full mass. Of course, evaporating the disk 
is a purely numerical contrivance, but this allows us to test for halo distribution function 
changes directly. Although classical mechanics are time- reversible, the random phases of any 
scattered orbits ensure that simply evaporating the central mass is not enough to return to 
the initial configuration. This would only be possible if we had a perfect integrator and if 
we had reversed all velocities, which we did not do. 



We formed prolate/triaxial halos via mergers, as described in iMoore et al.l (120041 ). The 



ini tial spherical halos were g enerated from a distribution function using the method described 



m 



Kazantzidis et al.l (120041 ) with the added refinement that each halo is composed of two 
mass species arranged on shells. The outer shell has more massive particles than the inner 
one, in order to increase the effective resolution in the central parts. Our model halo A was 
generated by the head-on merger of two prolate halos, themselves the product of a binary 
merger of spherical systems. The first merger placed the concentration c = 10 halos 800 
kpc apart approaching each other at 50 kms~^, while the second merger starts with the 
remnant at rest, 400 kpc from an identical copy. The resulting halo is highly prolate with a 
mild triaxiality. Halo model B was produced by the merger of two spherical halos starting at 
rest, 800 kpc apart. Both halos A and B consist of 4 x 10^ particles. The outer particles are 
~ 18 times more massive in halo A and ~ 5 times more massive in halo B. A large part of 
the segregation by particle mass per sists after the mergers and the small radius regions are 



dominated by low mass particles (cf. lDehnenll2005l ). Figure [T] shows the particle segregation 



in the case of halo A. We used a softening parameter e = 0.1 kpc for all halo particles, 
although we have verified that using a larger softening, e = 1 kpc, for the more massive 
species does not change our results. Our force resolution was chosen to be smaller than the 
vertical scale of the disk, thereby resolving short-range forces. 



Once we produced the prolate/triaxial halos, we inserted a disk of particles that remains 
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Fig. 1. — Spherically averaged density profile of dark matter particles in halo A before any 
baryons are introduced. The sohd hne is the full density profile, the dashed line is for the 
lower mass species, while the dotted line is for the higher mass species. 
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rigid throughout the experiments. In all mergers we have been careful to either give the halo 
no angular momentum, or to place the disk's symmetry axis along the angular monaentum 
of the halo since otherwise additional evolution would result (jPebattista fc Sellwoodlll999l ). 
The disks are composed of 3007^ equal-mass particles each with a softening e = 60 — 100 pc. 
The disk distribution was, in all cases, exponential with scale length and Gaussian scale- 
height Zd/Rd = 0.05. The disks were placed at various orientations within the halos. We 
refer to these experiments by the halo axis along which the disk's symmetry axis is aligned: 
in "short-axis" (S) experiments, the symmetry axis of the disk is parallel to the short axis 
of the halo, while in "long-axis" (L) experiments, the symmetry axis of the disk is along 
the halo's major axis. If the halo is triaxial, then an "intermediate-axis" (I) experiment has 
the disk minor axis parallel to the halo's intermediate axis. Initially, the disk has negligible 
mass, but this grows adiabatically linearly over time to a mass during a time tg. After 
this time, we slowly evaporated it during a time tg. Thus, 



0<t<t„ 



Mdisk{t) 



Mfe (1 - tg<t<tg + 



From t = Q to tg + tf. the halo particles are free to move and achieve equilibrium with the 
disk as its mass changes, but all disk particles are frozen in place. Since a triaxial global 
potential leads to elliptical disks forming, we include one simulation with an elliptical disk. 

Another key assumption in these simulations is that the disks form without much trans- 
fer of angular momentum to the halo. While formation of realistic galaxies requires that 
baryons conserve most of their angular momentum (e.g., ISommer-Larsen et al.lll999l). gas 
cond ensation onto subhalos results in angular momentum transfer to the halo (INavarro fc Steinmet2 
19971 ). We therefore present experiments in which a few softened particles were introduced, 
with a mass grown in the same way. We refer to these experiments by the label "P" sub- 
scripted by "f" for particles frozen in place and by "1" for live particles free to move. Lastly, 
we present one simulation, BAl, in which the disk at tg is replaced by live particles and 
evolved for 10 Gyr before evaporating the disk. We set initial disk particle velocities for 
a constant Toomre-Q = 1.5. In setting up the disk kinematics we azimuthally averaged 
radial and vertical forces. Thus, our disk is initially not in perfect equilibrium, but was close 
enough that it quickly settled to a new equilibrium. 

All the sirnulatio ns in this paper, which are listed in Table [H were evolved with PKD- 
GRAV (jStadeJ |200l|) , an efficient, multi-stepping, parallel treecode. 
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Run 


Halo 


^"200 


M200 


M, 


Rd 




te 






[kpc] 


[10^2 M^] 


[10^1 M^] 


[kpc] 


[Gyr] 


[Gyr] 


SAl 


A 


215 


4.5 


1.75 


3.0 


5 


2.5 


SA2 


A 


215 


4.5 


5.25 


3.0 


5 


2.5 


SA3 


A 


215 


4.5 


1.75 


1.5 


5 


2.5 


lAl 


A 


215 


4.5 


1.75 


3.0 


5 


2.5 


LAI 


A 


215 


4.5 


1.75 


3.0 


5 


2.5 


LBl 


B 


106 


0.65 


1.05 


3.0 


15 


7.0 


TAl 


A 


215 


4.5 


1.75 


3.0 


5 


2.5 


EAl 


A 


215 


4.5 


1.75 


3.0 


5 


2.5 


BAl 


A 


215 


4.5 


0.52 


3.0 


1.5 


2.5 


P^Al 


A 


215 


4.5 


1.75 


0.5 


5 


2.5 


PiA2 


A 


215 


4.5 


1.75 


5.0 


5 


2.5 


PiBl 


B 


106 


0.65 


0.7 


3.0 


10 


4 


P/B2 


B 


106 


0.65 


0.7 


3.0 


10 


4 


P^BS 


B 


106 


0.65 


0.35 


0.1 


5 


5 



Table 1: The simulations in this paper. For the particle simulations (PiAl-P;B3), refers 
to the softening of the particle(s). For runs P^Al and P/A2, the value of Mb refers to the 
combined mass of all the satellite particles at tg. The disk in run TAl is tilted by 30° 
relative to the one in SAl whereas the disk in EAl is elliptical. In run BAl, the disk at tg 
was replaced by live particles and evolved for 10 Gyr (during which time a bar formed and 
then was destroyed), before being evaporated. 



2.1. Measuring halo shapes 



To measure the axis ratios c/a and b/a we adopt a method based on iKata (Il99ll ) that 
uses the eigenvalues of the (unweighted) moment of inertia tensor /. For each bin of N 
particles we computed lij as follows: 



We then diagonalize I and calculate 

b/a= A/J22/X] 



11 



and c/a = a/Tss/Jh, 



(2) 



(3) 



where the Xjj's are the eigenvalues of I and Xn > X22 > X33. A useful pararneter f or expressing 
shape is the triaxiality parameter T = (a^ — 6^)/(a^ — c^) (jPranx et al. 1991 ). The cases 
T = and T = 1 correspond to oblate and prolate shapes, respectively, while T = 0.5 is the 
maximally triaxial case. 



We measured shapes in shells of fixed semi-major axis widths around the center of the 



syste m. Thus, these shape measuremen ts are differe ntial, rather than integrated (c/. iKatz 



199lh . We use the iterative procedure of iKata (jl99ll ) in which the convergence criterion is a 



variation in axis ratios by < 0.01%. In each iteration the semi major axis of the shell is held 
fixed; a particle is included in the calculation of lij if qio < q < qhi, where q is the ellipsoidal 
radius defined as 



We used shell widths qhi — qio = 5 kpc for all models. 

The center of the system is taken to be the center of mass of a sphere of radius 1 kpc 
centered on the minimum of the potential and is fixed for all shells. Tests performed in 
which the center of mass was allowed to vary by up to 0.5 kpc show less than 5% variation 
in the axis ratios past 10 kpc. Tests in which the limits of each shell were reduced by half 
around the average radius of that shell gave axis ratios that varied from the full resolution 
results by up to 0.08 in the worst cases and by less than 0.05 for most runs. Shells were not 
prevented from overlapping; as a result, some particles are sampled in more than one shell. 
We have verified that this does not bias our shape estimates through the experiments with 
the halved shell widths, where the shells never intersect. When the number of particles in 
the central shell is less than lOK, then convergence is not reached after 20 iterations, or the 
axis ratios varies by as much as 20%, so we take this number to be a reasonable cutoff for 
the reliability of these innermost shells and ignore shells with less particles. This occurs in 
only two cases, and in general most inner shells have > 25i^' particles, which we find to be 
more than enough to ensure consistent measurements with our method. 



3. Central Massive Objects 
3.1. Short- and Intermediate- Axis Experiments 



In run SAl we grew a disk inside halo A with the minor axes of the disk and halo 
aligned. This orientation is a natural one for disks to form in since simulatio ns have shown 
that the angular momenta of halos are aligned w ith their minor axes {e.g. IWarren et al. 
19921 : iPorciani et al.l 120021 : iFaltenbacher et al.ll2005l ). Once the disk is grown to its full mass 
it dominates the inner rotation curve (see Fig [2l left). The shape evolution of this highly 
prolate, mildly triaxial halo is shown in Figure (HH- The effect of the massive disk on the 
halo shape is large: having started out with b/a < 0.45 it becomes much rounder in the 
plane of the disk (b/a > 0.6 to 40 kpc, i.e. ~ 0.2 r2oo, the disk constituting 17% of the 
mass within this radius), as shown by the blue lines. The change in shape perpendicular to 



- 9 - 



the disk is more modest and the inner halo becomes significantly more triaxial than when 
it started out. In the disk plane, the combined potential starts out very elliptical, e$ ~ 0.4, 
and becomes quite round, with e$ < 0.1 over the entire extent of the disk (Fig. H]). This e$ 
i s sufficiently small to be consistent with the observed scatter in the TuUy-Fisher relation 



( iFranx fc de Zeeuwi Il992l ). even without the additional axisymmetrization of the potential 
that would be caused by the disk's orthogonal response. Once the disk is evaporated, the 
resulting halo shape, shown by the green lines in Figure [3^, is very similar to its original 
shape, with the net increase in both h/a and c/a being < 0.1 throughout the inner 100 kpc. 
The final triaxiality is barely changed from the starting one, despite the fact that the inner 
halo was almost maximally triaxial at tg. 

Likewise, the final density and anisotropy profiles, shown in Figure are not signifi- 
cantly changed, despite the factor of ~ 3.7 increase in halo central density at tg. The halo 
anisotropy, [3 = 1 — cr^^/cr^, where = |(cr| + cr^), starts out /? — 0, grows to P > 0.2 at tg 
{i.e. becomes radially anisotropic), and returns to /3 ~ at + tg- If box orbits had been 
destroyed to an y significant extent, we would have seen instead an increase in tangential 



anisotropy (e.g.. iHoUey-Bockelmann et al.ll2002l ) 



The small difference between t = and tg + in halo shape and kinematics suggests 
that the halo phase space distribution has not been grossly altered by the presence of the 
massive disk. There is little evidence for a substantial amount of box orbit scattering, and 
any chaos induced has to be quite mild. All this is true despite the quite large change in 
halo shape and kinematics when the disk is at full mass. 

Figure [3^ {yellow and red lines) also shows the evolution when we left the disk at full 
mass for 5 Gyr before evaporating it. The halo becomes slightly rounder (by about 6/a < 0.1) 
at all radii both after the additional 5 Gyr and once the disk is evaporated. The difference 
is largest inside ~ 30 kpc where the final b/a and c/a are about 0.1 larger at the end of 
the simulation. This difference must be due to orbit scattering (either physical or purely 
numerical); the fact that the difference between these two runs is so much smaller than that 
between t = and tg implies that scattering has only a mild effect on the halo shape. The 
global shape change at tg can therefore be attributed to orbit deformation. 

Model SAl had a disk with = 3 kpc a nd a baryon-to-dark matter fraction, = 0.039, 



consistent with estimates for local galaxies (j Jimenez et al.ll2003l ). A more massive or more 



compact galaxy may lead to greater scattering. We explored to what extent larger fh or 
smaller affect the halo shape in two further simulations. Run SA2 increased by a 
factor of 3 while keeping R^ fixed. The halo shape is changed significantly all the way out 
to r2oo once the disk is evaporated, but remains quite prolate, with b/a < 0.6 and c/a < 0.5, 
as can be seen in Figure [3]d. In contrast, at tg the halo has 0.5 < b/a < 1.0 within the 
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i nner 50 kpc. Even w ith this high ~ 0.12, or ~ 70% of the full cosmic baryon fraction 



(ISpergel fc et al.ll2006l ). the irreversible change to the halo shape is < 50% of the full change 
at tg out to 100 kpc. Run SA3 instead set = 1.5 kpc, keeping the ratio z^/Rd fixed (and 
decreasing all softenings appropriately). The evolution in this case is shown in Figure [3t; 
as in run SAl, although the halo at t = tg is substantially rounder than at the start, after 
the disk is evaporated the halo recovers most of its original shape. Of runs SA2 and SA3, 
making the disk more massive (SA2) produced a larger irreversible change in the halo than 
did making it more compact. 

Run lAl explored whether having the disk orthogonal to the intermediate axis makes a 
significant difference to the halo distribution, with all other parameters as in run SAl. The 
resulting shape evolution is presented in Figure [31i. In this simulation the halo at tg remained 
more elongated than that in run SAl despite having the same Mf,. The axis ratios of the 
halo at tg cross over at ~ 30 kpc, where the halo's flattening orthogonal to the disk causes 
the minor axis to switch from the disk plane to the orthogonal direction. Once the disk is 
evaporated, the halo ends very nearly axisymmetric in cross section in this inner region but 
continues to be highly prolate. As in run SAl, the net change in halo shape is relatively 
small at tg + te- 



3.2. Long-Axis Experiments 



In run LAl we placed the disk with its symmetry axis along the long axis of the halo. 
This orienta tion has been sugges ted to be favored by the distribution of satellit es around the 
Milky Wav (IZentner et al.l 120051 1 and by the Sagittarius dwarf tidal stream (IHelmil (120041 ) 
but see iFellhauer et al.l (120061 ) for a different view). Other than the disk's orientation, the 
parameters of this model are identical to those of run SAl. As in that model, the halo in 
run LAl is signiflcantly deformed to large radius by the growing disk, but it recovers its 
shape nearly completely once the disk is evaporated, as shown in Figure [3^. Likewise, the 
spherically averaged kinematic evolution of run LAl is indistinguishable from that of SAl, 
as seen in Figure [51 

A unique characteristic of the evolution in long-axis experiments is their tendency for 
the major axis of the inner halo to switch orientation by 90° into the disk plane once the 
disk grows sufficiently massive. For run LAl this is evident in Figure [3^, which shows that 
the halo is axisymmetric at ~ 20 kpc (the solid blue line approaches 6/a ^ 1, while the dot- 
dashed blue line approaches T ~ 0) but is quite prolate-triaxial at smaller radii. Major axis 
fiips are more clearly illustrated by the nearly axisymmetric halo B. In run LBl, as Mdisk 
increases, shells of the prolate inner halo become spherically symmetric. Further increase in 
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Mdisk then leads to the shell becoming not only flatter vertically (relative to the disk) but 
also acquiring an elongation with its major axis in the plane of the disk, i.e. the symmetry 
axis of the inner halo flips by 90° and becomes orthogonal to that of the outer halo (see Figs. 
[n]and[7j). The direction along which this reorientation occurs is not random since the halo is 
initially not perfectly axisymmetric on large scales. Continued increase in Mdisk causes the 
symmetry axis to flip orientation to larger radii, eventually saturating at ~ 10 kpc. The 
halo orientation flips do not occur when the disk is replaced by a point particle. The dotted 
gray line in Figure [7] shows the orientation of the major axis in run P;B1, at a time when its 
mass is the same as that in LBl; no flip in the major axis direction can be seen. 



3.3. Inclined Disk 



Sim ulations have found that t he angular morn enta of the halo and gas need not be 
aligned (Ivan den Bosch et al.l 120021 : IChen et al.ll2003r). although the inn er halo and the disk 
that would form settle to a common plane (IDubinski &: Kuijkeru Il995l ). Therefore we also 
explored the effect of a disk inclined relative to the main plane of the dark matter halo. 
Run TAl is based on run SAl but with the disk inclined by 30° about the ?/— axis (inter- 
mediate axis), all other parameters being the same. The resulting evolution is virtually 
indistinguishable from that in run SAl at all times, as shown in Figure [S^. 



3.4. Elliptical Disk 

A disk forming in an elliptic al potential becomes elo ngated with its major axis orthog- 



onal to that of the potential (e.g., iGerhard fc Vietrilll986l ). In run EAl we replaced the disk 



in run SAl by an elliptical disk, with its long axis along the ?/— axis. We obtain this oval 
disk by shrinking the x coordinates (parallel to the halo major axis) of all disk particles by 
a factor of 0.75; i.e., the ellipticity of the disk density was e = l — h/a = 0.25. The degree 
to which elliptical disks change halo distributions is over estimated by this simulation since 
the disk ellipticity is quite large for a massive disk. We then evolved this system identically 
to run SAl, keeping the disk flxed in place. Figure [Hb compares the shape evolution with 
that in run SAl. At tg + the halo is left significantly rounder within ~ 20 kpc than it was 
in run SAl, but beyond that the evolution is very similar. 



3.5. Central Softened Point Masses 



All of the experiments described above had rigid disks frozen in place. While we have 
been careful to recenter the halo in position and velocity after the mergers and before growing 
the disks, some residual relative motion of the inner and outer parts of the halos remained. 
This motion is damped as the mass of the disk increases, possibly causing some scattering of 
orbits. In order to test for artifacts associated with such damping, we resorted to simulations 
with only a single baryonic particle and compared the evolution when the particle is free to 
move (P^Bl) and when it is frozen in place (P/B2). The rotation curve at tg is shown in 
Figure [21 Figure [3k shows their shape evolution; in both cases, the shape is largely recovered 
at tg + te- If anything, P/Bl is very slightly rounder at 10 kpc than P/B2 compared with 
the initial halo. Figure [5] shows that their density and kinematic evolution also is largely 
reversible. Thus, our use of rigid disks nailed in place could not have induced much artificial 
orbit scattering. 



3.6. Ultrahard particle 

In runs P^Bl and P/B2 the growing particle had a softening of e = 3 kpc, a reasonable 
size for a galaxy. In run P;B3 we decreased the softening length of the particle to 100 pc and 
Mb by half. Despite the smaller the final halo after tg + te remains substantially rounder 
inside 20 kpc than in those runs (but is largely recovered at larger radii). If the central 
particle were a black hole, its sphere of influence assuming = 100 kms~^ from t = would 
be ~ 15 kpc. This is comparable to the radius out to which the particle irreversibly alters 
the shape of the halo. Likewise, the halo mass within 20 kpc is comparable to that of the 
central particle: at t = 0, the halo mass within this radius is 4M;,. Unlike runs P^Bl and 
PjB2, Figure [3 shows that the kinematic evolution is not reversible, and the halo of run 
P;B3 becomes significantly tangentially anisotropic, as expected if box orbits are destroyed. 
Despite the different final state. Figure [3)i also shows that the halo shape at tg is not much 
different from that in P;B1 and P/B2, implying that halo shape change is not dominated by 
scattering. Whereas run SA3 with = 1.5 kpc, which is not unreasonably small for most 
galaxies, did not significantly cause box orbit destruction, the ~ 10 times more centrally 
concentrated run PiB3 is able to cause a large irreversible change to the halo shape out to 
~ 0.3 r2oo- However, Rd/r2oo — 0.001 is unrealistically small. 
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4. Angular Momentum Transport 

We next explore the effect of angular momentum transfer to the halo. Such transfer 
is irreversible so the change inflicted on the halo must also be irreversible. How strongly 
the halo distribution is changed depends on the mechanism by which angular momentum is 
transferred. If via bars or spirals then we may expect that the changes are mostly at small 
radius. If angular momentum is transferred by satellite galaxies, however, then the effect on 
the halo is likely to be much more widespread. 



4.1. Live Barred Disk 



In run BAl we evolved a model with a live disk for 10 Gyr after tg before evaporating 
the disk. The initial system was similar to run SAl but with only 30% of its Mf,. We chose 
this lower mass because the same mass as SAl leads to a long-lived bar, whereas we are 
interested in forming a bar that gets destroyed in ord er to be able to evapo rate the disk. A 
bar quickly formed and was subsequently destroyed (IBerentzen et al.ll2006l ). After 10 Gyr, 
we fixed the disk particles in place and evaporated the disk. Very little of the inner halo 
shape is recovered after the disk is evaporated. In the inner ~ 20 kpc the halo remains 
rounder than at t = 0, with both b/a and c/a larger by > 0.1. The irreversib le change in the 
halo is associated with the tr ansfer of angular momentum from disk to halo (jWeinberglll985l : 
Debattista &: Sellwoodlll998l ). Run BAl produces a comparable change in the inner 20 kpc 
of the halo as did the 10 times more massive run SA2, but the shape change is mu ch smaller 
farth er out. Since the bar transfers angular momentum to the halo at resonances (jWeinberg 
19851 ). and the strongest of these are at smaller radii (most of the angular momentum gained 
by the halo is within the inner ~ 10 kpc), this accounts for the relatively small radial extent 
of the halo shape change. 



4.2. Satellites 



Baryons need not cool directly onto the central disk but onto satellites inste ad. The pres 
ence of large numbers of dark satellites is one of the main predictions of CDM (IMoore et al. 
19991 : iKlypin et al.l 1 19991 : iGhigna et al.l |2000| ) . As they sink, satellites lose angular momen- 
tum to the halo; in the process box orbits may be scattered. We explored this evolution 
with models P;A1 and PiA2. Starting with halo A, we selected 10 particles that stay within 
200 kpc but otherwise at random, and adiabatically increased each of their masses to give 
the same total baryonic mass as in run SAl. We grew these satellites to full mass and then 
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evaporated them. Each satelhte had e = 0.5 kpc in P;A1 and 5 kpc in P;A2. Since we only 
used softened point particles as satellites, which cannot be tidally stripped, their effect on 
the halo is larger than it would be in nature. Of the 10 satellites, only one remained at 
r > 50 kpc, the rest having fallen to i? < 25 kpc by the end of the simulation. The evolution 
of the halo shape in these two models is presented in Figure [3]F. After the particles reach 
their full mass, the halo of run P/Al is about as round as that in run SAl. However, the 
halo does not recover much of its original shape after the particles are evaporated. Clearly, 
the distribution function of the halo has been altered to a large extent. Figure shows 
that angular momentum tra nsferred by baryons to the halo can erase the cusp, in agree- 



ment with previous res ults (iTonini et al.l l2006l : iMashchenko et al.l l2006l : iRead et al.l 12006 



Weinberg fc Katzll2007l ) although the contraction caused by the growing central mass masks 



the core. The halo shape change at tg + tf, is very similar in the two runs because of quite 
similar angular momentum absorbed by the halo. Because of the difference in softenings, 
the different baryonic potential at the halo center at tg accounts for the b/a ^ 0.1 difference 
in shape, with the softer potential supporting the more elongated shape. 



5. Orbital Evolution 



We explored directly the evolution of the orbital character of the models by considering 
a subsample of 1000 particles in run SAl and following their orbits at various points in the 
simulation. The 1000 particles were randomly chosen from the t = distribution such that 
they were inside r = 200 kpc. We then integrated their motion as test particles while holding 
all the other particles fixed in place. We used a fixed timestep of 0.1 Myr and integrated 
for 15 Gyr, storing the phase space coordinates of each test particle every 1 Myr. For the 
same 1000 particles, we carried out this operation at t = 0, tg, and tg + te- Because we froze 
the background potential, in effect we have computed the orbital character of the particles 
at these three times. The fact that we integrated for 15 Gyr ensures that we have sufficient 
points on each orbit to properly characterize it. In this paper we demonstrate with a few 
examples that a large fraction of box-like orbits at t = return to very similar box-like orbits 
at tg + te, showing that deformation, not transformation, is responsible for shape change in 
most cases. We do this by presenting their configuration space projection at each of the 
three different times. Such an analysis cannot distinguish between box orbits and mildly 
chaotic, elongated orbits but this is unimportant anyway for our present purposes since we 
have integrated for over a Hubble time. If they are mildly chaotic, they can still support 
a triaxial halo. A full analysis of the orbital structure using more sophisticated techniques 
will be presented elsewhere. 
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Of the 1000 orbits, we start by presenting nine particles, initially on boxlike orbits, 
defined such that, at t = 0, they (1) remain inside 25 kpc, (2) do not have a fixed sense 
of rotation relative to any of the three major axes, (3) reach a radius of at least 10 kpc, 
and (4) get within 0.2 kpc of the center. The evolution of many of the other 991 orbits is 
qualitatively similar to that of the nine presented here. Figure [9] projects these orbits onto 
the halo symmetry planes, where the x-axis is the halo's major axis and the ^-axis is the 
disk's symmetry axis. Most orbits at tg + tg are quite similar to what they looked like at 
t = 0. None of the orbits seem strongly chaotic, neither at tg nor at tg + tg, although they 
may be weakly chaotic. Moreover, most orbits retain a box-like shape at tg, but have a 
significantly rounder shape than those at t = 0. At tg, three of the initially box-like orbits 
become round (orbits "a", "f and "h"); of these only orbit h changes character completely, 
becoming a loop orbit. Some of the orbits have a slight banana shape; in the full sample 
of orbits we found many cases of strongly banana-shaped orbits. These had a tendency to 
become more planar but are still distinctly elongated at tg + te- In a few cases we also found 
the opposite occurring — slightly banana orbits becoming more strongly curved — but this 
was less common. Of the box-like orbits in Figure [9] some are rounder in the {x, y) plane at 
tg + te {e.g., orbits a, f and i), but some are rounder at t = {e.g., orbits c and e), suggesting 
that differences in shape are due to scattering. What little difference in orbit shape occurs 
between t = and tg + tg can probably be attributed to numerical noise. Most importantly, 
while there is a clear orbit shape deformation at tg, little orbital transformation has occurred. 

We quantify the orbital deformation of the sample of 1000 particles by plotting in 
Figure [TUl(T,,/cr^.. where = J2t-^t similarly for and the sum is over timesteps. The 
significantly rounder shape of orbits at tg than t = is apparent, with the vast majority of 
orbits initially aligned with the halo having (Ty/ax closer to unity at the later time. Orbits 
initially elongated along the halo's minor axis, as well as orbits initially rounder than ay/a^ ^ 
0.6 end up round, with ay/a^ — 1. Instead at tg + te the orbits tend to return to their initial 
elongation, especially for the most elongated orbits. The right panel shows the distribution 
of Oyjox, orbits become substantially rounder at tg but the population as a whole recovers 
the original distribution to a large extent once the disk is evaporated. 

A comparison of orbital evolution in runs P/B2 and P^BS provides an example of orbit 
transformation. We again selected 1000 orbits from particles within the inner 200 kpc of halo 
B at t = 0. We integrated their orbits as above but we used the smaller timestep 6t = 10^ 
years in the case of P^BS; for P/B2 we use tg/2 when the central particle has the same mass as 
at tg in model P/BS. Although we are comparing the two models at the same central particle 
mass, the orbits at time tg + tg in model P/B2 were computed after the central particle was 
evaporated from a mass twice that reached in P^BS. As before, we present in Figure [TT] nine 
orbits that at t = are box-like. These boxlike orbits in P;B3 are more often transformed 
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than those in run P/B2. In model P/B2, only one orbit (orbit f) appears to have changed 
substantially at the end of the simulation, while orbits a and d become fish orbits (although 
they may have been librating about fish orbits at t = 0). Orbit b is largely unchanged and 
the remaining orbits are all boxlike. In model P;B3, orbit e is changed about as much as 
orbit f in P/B2. However, four of the nine orbits, f-i, are very strongly transformed by tg + te 
and are no longer able to support a triaxial shape. Figure [12] compares the distribution of 
all 1000 particles and a depletion of elongated orbits is evident in P^BS compared with run 
PfB2. 



6. Discussion 
6.1. Timesteps 

We also performed a number of tests of the numerics to veri f y that our results are 



robust. The main concern is the timestep used. IShen fc SellwoodI (120041 ) found that too 
large timesteps result in bars being destroyed too easily, because orbits are not followed 
accurately near the central mass. Our simulations used multi-stepping. With a base timestep 
At, particles move on timesteps At/2"', where n is the rung level satisfying the condition 
dt = At/2" < r](e/aY^'^, and e is the particle's softening, a is its acceleration and r] a 
tolerance parameter. We used 1] = 0.2, a conservative value; with a base timestep At = 5 
Myr, simulation SAl at tg had a range of timesteps down to 5/2^ = 0.16 Myr. If instead 
we set rj = 2 the timestep distribution only reaches to 5/2 = 2.5 Myr. The effect of these 
larger timesteps, shown in Figure [T31 is manifest at r < 20 kpc, which remains significantly 
rounder at tg + te than when rj = 0.2. The quite modest net shape change in our simulations 
implies that the timesteps we used were sufficiently small to correctly follow the evolution 
near the center. 



6.2. Evidence against box orbit destruction 



We have demonstrated that the substantial axisymmetrization caused by a disk growing 
inside a dark matter halo is largely, although n ot wholly, reversible . If triaxial halos become 



rounder because of box orbit destruction {e.g. iMaccio et al.l 120071 ) . then these orbits would 
have to be repopulated in order for the halo to recover its original shape once the disk is 
evaporated. Apart from being unlikely for such a highly ordered system as a triaxial halo, this 
interpretation is not supported by our orbital analysis, which shows that box-like orbits are 
deformed by realistic disks but not transformed by scattering into new orbits. The strongest 
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evidence that deformation is a more important process than transformation comes from the 
abihty of particles to return, after the disk is evaporated, to nearly the same orbits as they 
started from. This makes it implausible that the shape change is due to a large increase 
in strong chaos. Instead we find that realistic axisymmetric di sks are not concentrated 



or ma ssive enough to cause substantial chaos. In a similar vein, iHoUey-Bockelmann et al. 



( I2OO2I ) found that chaos in triaxial ellipticals is induced by the black hole only, not by the 
stellar cusp surrounding it. This need not mean, however, that chaos is not enhanced by the 
disk. It could well be that orbits are becoming weakly chaotic but do not diffuse sufficiently 
on a Hubble timescale to significantly weaken triaxiality. The shape evolution due to baryons 
cooling onto a central galaxy can therefore be computed from adiabatic invariants. 

The primary role of orbit deformation over transformation is also indicated by the 
much smaller effects of scattering when it is clear that scattering has occurred. In run 
SAl scattering occurred when the disk was maintained at full mass for an additional 5 Gyr 
before being evaporated. This only caused a small additional change in the final halo shape. 
Furthermore, comparing runs P;B1 and P/BS, we find that the halo shape at tg is rather 
similar, despite the fact that in run F^BS box orbits are significantly destroyed. The same is 
true for runs SAl and F^Al. These examples directly illustrate that box orbit destruction 
is a much smaller factor in halo shape changes than is orbit deformation. 



Weak chaos is consistent with the orbital characterization of iMaccio et al.l (120071 ) but 
does not support their claim that the shape change is largely due to enhanced chaos from 
the central baryonic mass. Another possibility is that in their simulations chaos could have 
been caused by gas cooling inside subhalos, which we showed leads to a substantial change 
in the phase space distribution of the halo. 



6.3. Implications for galaxy formation simulations 



If baryons cool onto substructures within the halo then box orbits are very efficiently 
destroyed. Since the evolution of disks can be strongly infiuenced by halo triaxiality {e.g. 
Ideta fc Hozumil I2OOOI ) . any process that artificially reduces tri axiality can lead to biases in 
the properties of galaxies forming in cosmological simulations. lAgertz et al.l (120071 ) demon- 
strate that the evolution of gas blobs in smoothed particle hydrodynamics (SFH) simulations 
is different from that found in Eulerian gas codes. They interpreted this difference as be- 
ing due to the unphysically poor mixing of traditional SFH, which allows blobs to survive 
longer in SFH. Moreover, satellites in cosmological simulations tend to have denser, more 
concentrated gas component s than their real co unterparts, which makes them harder to strip 
by ram pressure and tides (IMayer et al.l 120071 ) and thus more likely to artificially enhance 
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halo shape changes. Baryonic cooling inside substructures is strong ly suppressed by feed- 
back from supernovae (jPekel &: Silklll986l : iGovernato et al.ll2004l . 120071 ) . If simulations do not 
treat feedback properly or have low resolution, then baryons may condense into concentrated 
substructures biasing the global evolution of simulated galaxies. 



6.4. Summary 



Our results can be summarized as follows: 



The adiabatic growth of disks with realistic sizes and masses inside prolate/triaxial 
halos leads to a large change in the shape of the halo. Axis ratios change by > 0.2 out 
to roughly 0.5 r2oo- The growth of the disk drives the halo kinematics to larger radial 
anisotropy. The midplane global potential ell ipticity is less than 0.1, consistent with 
the small scatter in the TuUy-Fisher relation ( Franx fc de Zeeuwlll992l ). 



2. Despite these large changes, the underlying phase space distribution is not grossly 
altered, as we verified by artificially evaporating the disk and recovering, to a large ex- 
tent, the original halo. The irreversible change in final halo structure is larger for more 
massive or more centrally concentrated disks, but is still a relatively small fraction of 
the total shape change when the dis k is at full mass. As in the case of black holes at the 



centers of cuspy elliptical galaxies (iGerhard fc Binneylll985t iHoUey-Bockelmann et al. 



2OO2I ). the bulk of the irreversible halo shape change occurs in the inner region of the 



galaxy. This small irreversible shape change is driven by orbit scattering. 

3. Box orbit destruction cannot be the right interpretation for the shape change caused 
by disk growth. Such a process is not reversible but we found that a large fraction of 
particles on box-like orbits individually return to very similar orbits after the disk is 
evaporated. At most only mild chaos is induced. Instead we find that box orbits be- 
come deformed by the growing disk, but retain their character, and this seems sufficient 
to explain the change in shape. In the absence of angular momentum transport or ex- 
treme mass/concentration galaxies, very little of the quite large shape change that dark 
matter halos undergo as baryons condense inside them is due to box orbit destruction. 
As a result, shape change can be well approximated by adiabatic invariants. 

4. Very concentrated structures do lead to scattering and to large irreversible changes in 
halo shape and kinematics; we found that scales of ~ 100 pc are needed to accomplish 
this. The irreversible shape change was then restricted to the sphere of influence of 
this pointlike mass: r ~ GM/a"^. Massive disks are also able to change the halo 
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structure irreversibly, but the mass required, ~ 70% of the cosmic baryon fraction, is 
quite high. Even in such cases of scattering, the halo shape with the baryons at full 
mass is not much different from similar simulations with little scattering, suggesting 
that it is orbital deformation in the first place that changes the shape of the halo. 

If baryons transport angular momentum to the halo, a large irreversible change in 
halo shape and kinematics occurs. Such transfers can occur either because of gas 
condensing in satellites or nonaxisymmetric structures forming in the disk. Even quite 
low mass disks are able to alter the inner halo distribution if they can transport angular 
momentum to the halo. The effect of satellites can be art ificially large in sim ulations 
of galaxy formation because of the poor mi xing in SPH (lAgertz et al.l 120071 ) and the 



too highly concentrated satellites that form (IMayer et al. 



2003). 



6. When the disk minor axis and halo major axis are aligned, growth of the disk leads 
to an elongation within the plane of the disk, even when the initial halo is very nearly 
axisymmetric. 



This paper is based in part on work supported by the National Science Foundation un- 
der the following NSF programs: Partnerships for Advanced Computational Infrastructure, 
Distributed Terascale Facility (DTF), and Terascale Extensions: Enhancements to the Ex- 
tensible Terascale Facility. Many additional simulations were carried out at the University 
of Ziirich on zBox and at the Arctic Region Supercomputing Center. V. P. D. thanks the 
University of Ziirich for hospitality during part of this project. V. P. D. is supported by a 
Brooks Prize Fellowship in Astrophysics at the University of Washington and receives partial 
support from NSF ITR grant PHY-0205413. S. K. acknowledges support by the US Depart- 
ment of Energy through a KIPAC Fellowship at Stanford University and the Stanford Linear 
Accelerator Center. We thank Andrea Maccio and loannis Sideris for fruitful discussions and 
the anonymous referee for comments that helped to improve this paper. V. P. D. thanks 
Monica Valluri for discussions that helped improve this paper. 



REFERENCES 

Agertz, O., et al. 2007, MNRAS, 380, 963 



Allgood, B., Flores, R. A., Primack, J. R., Kravtsov, A. V., Wechsler, R. H., Faltenbacher, 
A., & Bullock, J. S. 2006, MNRAS, 367, 1781 



-20- 

Andersen, D. R., Bershady, M. A., Sparke, L. S., Gallagher, J. S., & Wilcots, E. M. 2001, 
ApJ, 551, L131 

Athanassoula, E., Lambert, J. C, & Dehnen, W. 2005, MNRAS, 363, 496 
Bailin, J., & Steinmetz, M. 2005, ApJ, 627, 647 

Bardeen, J. M., Bond, J. R., Kaiser, N., & Szalay, A. S. 1986, ApJ, 304, 15 

Barnes, E. L, & Sellwood, J. A. 2003, AJ, 125, 1164 

Barnes, J., & Efstathiou, G. 1987, ApJ, 319, 575 

Bartelmann, M., Steinmetz, M., & Weiss, A. 1995, A&A, 297, 1 

Berentzen, I., Shlosman, I., & Jogee, S. 2006, ApJ, 637, 582 

Buotc, D. A., & Canizares, C. R. 1994, ApJ, 427, 86 

Buote, D. A., Jeltema, T. E., Canizares, C. R., & Garmire, G. R 2002, ApJ, 577, 183 
Chen, D. N., Jing, Y. R, & Yoshikaw, K. 2003, ApJ, 597, 35 
Contopoulos, G. 1980, A&A, 81, 198 
Debattista, V. P. 2003, MNRAS, 342, 1194 

Debattista, V. R, Mayer, L., CaroUo, C. M., Moore, B., Wadsley, J., & Quinn, T. 2006, 
ApJ, 645, 209 

Debattista, V. P., & Sellwood, J. A. 1998, ApJ, 493, L5 
Debattista, V. P., & Sellwood, J. A. 1999, ApJ, 513, L107 
Dehnen, W. 2005, MNRAS, 360, 892 
Dekel, A., & Silk, J. 1986, ApJ, 303, 39 
Diehl, S., & Statler, T. S. 2007, ApJ, 668, 150 
Dubinski, J. 1994, ApJ, 431, 617 
Dubinski, J., & Carlberg, R. G. 1991, ApJ, 378, 496 
Dubinski, J., & Kuijken, K. 1995, ApJ, 442, 492 



Faltenbacher, A., AUgood, B., Gottlober, S., Yepes, G., & Hoffman, Y. 2005, MNRAS, 362, 
1099 

Fellfiauer, M., et al. 2006, ApJ, 651, 167 

Franx, M., & de Zeeuw, T. 1992, ApJ, 392, L47 

Franx, M., Illingwortfi, G., & de Zeeuw, T. 1991, ApJ, 383, 112 

Franx, M., van Gorkom, J. H., & de Zeeuw, T. 1994, ApJ, 436, 642 

Frenk, C. S., Wfiite, S. D. M., Davis, M., & Efstatfiiou, G. 1988, ApJ, 327, 507 

Fridman, T., & Merritt, D. 1997, AJ, 114, 1479 

Gerliard, O. E., & Binney, J. 1985, MNRAS, 216, 467 

Gerhard, O. E., & Vietri, M. 1986, MNRAS, 223, 377 

Gfiigna, S., Moore, B., Governato, F., Lake, G., Quinn, T., & Stadel, J. 2000, ApJ, 544, 616 
Governato, F., et al. 2004, ApJ, 607, 688 

Governato, F., Willman, B., Mayer, L., Brooks, A., Stinson, G., Valenzuela, O., Wadsley, 
J., & Quinn, T. 2007, MNRAS, 374, 1479 

Helmi, A. 2004, MNRAS, 351, 643 

Holley-Bockelmann, K., Mihos, J. C, Sigurdsson, S., Hernquist, L., & Norman, C. 2002, 
ApJ, 567, 817 

Huizinga, J. E., & van Albada, T. S. 1992, MNRAS, 254, 677 

Ibata, R., Lewis, G. F., Irwin, M., Totten, E., & Quinn, T. 2001, ApJ, 551, 294 

Ideta, M., & Hozumi, S. 2000, ApJ, 535, L91 

lodice, E., Arnaboldi, M., Bournaud, F., Combes, F., Sparke, L. S., van Driel, W., & 
Capaccioli, M. 2003, ApJ, 585, 730 

Jimenez, R., Verde, L., & Oli, S. R 2003, MNRAS, 339, 243 

Jing, Y. R, & Suto, Y. 2002, ApJ, 574, 538 

Jofmston, K. V., Law, D. R., & Majewski, S. R. 2005, ApJ, 619, 800 



-22- 

Jungman, G., Kamionkowski, M., & Griest, K. 1996, Phys. Rep., 267, 195 
Kalapotharakos, C, Voglis, N., & Contopoulos, G. 2004, A&A, 428, 905 
Katz, N. 1991, ApJ, 368, 325 

Kazantzidis, S., Kravtsov, A. V., Zentner, A. R., AUgood, B., Nagai, D., & Moore, B. 2004, 
ApJ, 611, L73 

Kazantzidis, S., Magorrian, J., & Moore, B. 2004, ApJ, 601, 37 

Klypin, A., Gottlober, S., Kravtsov, A. V., & Khokhlov, A. M. 1999, ApJ, 516, 530 

Kochanek, C. S. 1995, ApJ, 445, 559 

Koopmans, L. V. E., de Bruyn, A. G., & Jackson, N. 1998, MNRAS, 295, 534 
Kuijken, K., & Tremaine, S. 1994, ApJ, 421, 178 
Lake, G., & Norman, C. 1983, ApJ, 270, 51 

Maccio, A. V., Sideris, I., Miranda, M., Moore, B., & Jesseit, R. 2007, ArXiv e-prints, 704 
Mashchenko, S., Couchman, H. M. R, & Wadsley, J. 2006, Nature, 442, 539 
Mayer, L., Kazantzidis, S., Mastropietro, C., & Wadsley, J. 2007, Nature, 445, 738 
Merritt, D., & Poon, M. Y. 2004, ApJ, 606, 788 
Merritt, D., & Quinlan, G. D. 1998, ApJ, 498, 625 

Moore, B., Ghigna, S., Governato, F., Lake, G., Quinn, T., Stadel, J., & Tozzi, P. 1999, 
ApJ, 524, L19 

Moore, B., Kazantzidis, S., Diemand, J., & Stadel, J. 2004, MNRAS, 354, 522 

Navarro, J. F., & Steinmetz, M. 1997, ApJ, 478, 13 

Norman, C. A., May, A., & van Albada, T. S. 1985, ApJ, 296, 20 

Oguri, M., Lee, J., & Suto, Y. 2003, ApJ, 599, 7 

Poon, M. Y., & Merritt, D. 2002, ApJ, 568, L89 

Porciani, C., Dekel, A., & Hoffman, Y. 2002, MNRAS, 332, 325 



-23- 

Read, J. I., Goerdt, T., Moore, B., Pontzen, A. P., Stadel, J., & Lake, G. 2006, MNRAS, 
373, 1451 

Sackett, P. D., & Sparke, L. S. 1990, ApJ, 361, 408 

Schoenmakers, R. H. M., Franx, M., & de Zeeuw, P. T. 1997, MNRAS, 292, 349 
Schwarzschild, M. 1979, ApJ, 232, 236 

Schweizer, F., Whitmore, B. C, & Rubin, V. C. 1983, AJ, 88, 909 

Shen, J., & Sellwood, J. A. 2004, ApJ, 604, 614 

Sommer-Larsen, J., Gelato, S., & Vedel, H. 1999, ApJ, 519, 501 

Spergel, D. N., & et al. 2006, astro-ph/0603449 

Stadel, J. G. 2001, Ph.D. Thesis, University of Washington 

Statler, T. S. 1987, ApJ, 321, 113 

Tonini, C, Lapi, A., & Salucci, P. 2006, ApJ, 649, 591 

Udry, S., & Martinet, L. 1994, A&A, 281, 314 

Valluri, M., & Merritt, D. 1998, ApJ, 506, 686 

van den Bosch, F. C., Abel, T., Croft, R. A. C., Hernquist, L., & White, S. D. M. 2002, 
ApJ, 576, 21 

Warren, M. S., Quinn, P. J., Salmon, J. K., & Zurek, W. H. 1992, ApJ, 399, 405 

Weijmans, A.-M., Krajnovic, D., van de Ven, G., Oosterloo, T. A., Morganti, R., & de 
Zeeuw, P. T. 2007, ArXiv sprints, 711 

Weinberg, M. D. 1985, MNRAS, 213, 451 

Weinberg, M. D., & Katz, N. 2007, MNRAS, 375, 425 

Zentner, A. R., Kravtsov, A. V., Gnedin, O. Y., & Klypin, A. A. 2005, ApJ, 629, 219 



This preprint was prepared with the A AS macros v5.2. 



-24- 




5 10 15 20 5 10 15 20 



R [kpc] R [kpc] 

Fig. 2. — Azimuthally averaged rotation curves of models SAl (left) and PiBl/PjB2 {right) 
measured in tiie midplane. In both panels the solid line is the full rotation curve, the dashed 
line the contribution of the baryons, and the dot-dashed line the contribution from the halo. 
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Fig. 3.— Shape evolution in runs (a) SAl, {h) SA2, (c) SA3, (c?) lAl, (e) LAI, (/) P^Al and 
P/A2, i^g) P^Bl and P/B2 and (/i) P^BS. The sohd hues show 6/a, the dashed hues show 
c/a, and the dot-dashed hues show T (with scale indicated on the right-hand side of each 
panel). The black, blue, and green lines are at t = 0, tg, and tg -|- tg) respectively. In panel 
(a) black/yellow/red shows the evolution if, after t^, the disk is held at full mass for a further 
5 Gyr before it is evaporated. In panel (/), the standard colors are for P/Al, while PzA2 
is indicated in black/yellow/red. Likewise, in panel (a), the standard colors are for PfB2, 
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Fig. 4. — Ellipticity of the potential, e$, in the disk mid-plane of run SAl. The solid points 
are at t = 0, and the open points are at tg, with squares for halo only and circles for 
disk+halo. 
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Fig. 5. — Evolution of the spherically averaged density and kinematics in runs SAl [top left), 
P/Al [top right), P^Bl and P/B2 [bottom left), and PzB3 [bottom right). The black, blue, 
and green lines correspond to t = 0, tg, and tg + te- The red lines show tg + te for P/B2. The 
top panels shows the densities. In the middle panels the solid lines indicate ar, the dashed 
lines a^, and the dot-dashed lines a^. Here the z-axis from which the angle 9 is measured is 
the long axis of the halo. The bottom panels show the anisotropy parameter, /3. 
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Fig. 6. — Evolution of the inner halo in run LBl, seen in cross section in the disk plane, 
with only the region \z\ < 5 kpc shown, where the 2;-axis is the symmetry axis of both disk 
and halo. The panels show t = (/e/t), tg {middle) and tg + 1^. {right). The halo is initially 
axisymmetric, becomes elongated orthogonal to the symmetry axis at small radii at tg, and 
largely recovers its axisymmetry at tg + ig- 
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Fig. 7. — Major axis orientation evolution in run LBl. The solid black and gray lines are 
at t = and tg/3, the dashed black and gray lines are at 2tg/3 and tg and the black dotted 
line is at tg + tg. The gray dotted line shows run P;B1 at tg. 
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Fig. 8. — Shape evolution in runs (a) TAl, (6) EAl, and (c) BAl. The sohd hues show 6/a, 
the dashed hues show c/a and the dot-dashed hues show T (with scale indicated on the right 
hand side of each panel). The black, yellow, and red lines show t = 0, tg, and tg + 1^. For 
comparison, the equivalent evolution for run SAl is shown by the blue and green lines. 
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Fig. 9. — Sample of initially boxlike orbits in run SAl. In each set of panels, the top row is 
at i = 0, the middle row is at t^, and the bottom row is at tg + te] from left to right the panels 
show projections onto the {x, y), {x, z), and {y, z) planes. Each orbit has been integrated for 
15 Gyr from each of t = 0, tg and tg + te- 
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Fig. 10. — Deformation of orbits in run SAl. In the left panels, the dashed line shows the 
diagonal. The solid points show the nine particles of Fig. [91 In the right panel, the solid, 
dashed, and dot-dashed lines show the distributions at t = 0, tg, and tg + te, respectively. 
The thick lines are for the full distribution of 1000 particles, while the thin lines are for those 
particles that remain inside r = 50 kpc during the 15 Gyr integration starting at t = 0. 
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Fig. 11. — Sample of initially boxlike orbits in run P/B2 {left three columns) and P/B3 {right 
three columns). We show the evolution of the same nine starting orbits in the two different 
models. In each set of panels, the top row is at i = 0, the middle row is at tg/2 (P/B2) or tg 
(P;B3), and the bottom row is at tg + te] from left to right the panels show projections onto 
the {x,y), {x,z), and {y,z) planes. Each orbit has been integrated for 15 Gyr from each 
starting point. 
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Fig. 12.— Identical to Fig. [lObut for runs P/B2 (left) and P^BS (right). The final distri- 
bution of orbits is significantly depleted of elongated orbits in run P^BS compared with run 
P/B2. All nine orbits of Fig. [H] are above the diagonal in run P/B3 at tg + te but scatter 
about the diagonal in run P/B2. 
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Fig. 13. — Effect of timestep size on the evolution of run SAl. The black lines shows the 
shape at t = while the gray lines are for tg + t^, with the thick line the standard result 
with rj = 0.2 and the thin lines for rj = 2. As in Fig. [3l the solid lines show b/a while the 
dashed lines show c/a. The larger timesteps lead to an evolution that is less reversible. 



